Full-length transcriptome analysis of Ophioglossum vulgatum: effects of experimentally identified chloroplast gene clusters on expression and evolutionary patterns.

Genes with similar or related functions in chloroplasts are often arranged in close proximity, forming clusters on chromosomes. These clusters are transcribed coordinated to facilitate the expression of genes with specific function. Our previous study revealed a significant negative correlation between the chloroplast gene expression level of the rare medicinal fern Ophioglossum vulgatum and its evolutionary rates as well as selection pressure. Therefore, in this study, we employed a combination of SMRT and Illumina sequencing technology to analyze the full-length transcriptome sequencing of O. vulgatum for the first time. In particular, we experimentally identified gene clusters based on transcriptome data and investigated the effects of chloroplast gene clustering on expression and evolutionary patterns. The results revealed that the total sequenced data volume of the full-length transcriptome of O. vulgatum amounted to 71,950,652,163 bp, and 110 chloroplast genes received transcript coverage. Nine different types of gene clusters were experimentally identified in their transcripts. The chloroplast cluster genes may cause a decrease in non-synonymous substitution rate and selection pressure, as well as a reduction in transversion rate, transition rate, and their ratio. While expression levels of chloroplast cluster genes in leaf, sporangium, and stem would be relatively elevated. The Mann-Whitney U test indicated statistically significant in the selection pressure, sporangia and leaves groups (P < 0.05). We have contributed novel full-length transcriptome data resources for ferns, presenting new evidence on the effects of chloroplast gene clustering on expression land evolutionary patterns, and offering new theoretical support for transgenic research through gene clustering.

Chloroplast gene expression: Recent advances and perspectives.

Chloroplasts evolved from an ancient cyanobacterial endosymbiont more than 1.5 billion years ago. During subsequent coevolution with the nuclear genome, the chloroplast genome has remained independent, albeit strongly reduced, with its own transcriptional machinery and distinct features, such as chloroplast-specific innovations in gene expression and complicated post-transcriptional processing. Light activates the expression of chloroplast genes via mechanisms that optimize photosynthesis, minimize photodamage, and prioritize energy investments. Over the past few years, studies have moved from describing phases of chloroplast gene expression to exploring the underlying mechanisms. In this review, we focus on recent advances and emerging principles that govern chloroplast gene expression in land plants. We discuss engineering of pentatricopeptide repeat proteins and its biotechnological effects on chloroplast RNA research; new techniques for characterizing the molecular mechanisms of chloroplast gene expression; and important aspects of chloroplast gene expression for improving crop yield and stress tolerance. We also discuss biological and mechanistic questions that remain to be answered in the future.

Environmental and nuclear influences on microalgal chloroplast gene expression.

Microalgal chloroplasts, such as those of the model organism Chlamydomonas reinhardtii, are emerging as a new platform to produce recombinant proteins, including industrial enzymes, diagnostics, as well as animal and human therapeutics. Improving transgene expression and final recombinant protein yields, at laboratory and industrial scales, require optimization of both environmental and cellular factors. Most studies on C. reinhardtii have focused on optimization of cellular factors. Here, we review the regulatory influences of environmental factors, including light (cycle time, intensity, and quality), carbon source (CO2 and organic), and temperature. In particular, we summarize their influence via the redox state, cis-elements, and trans-factors on biomass and recombinant protein production to support the advancement of emerging large-scale light-driven biotechnology applications.

Transgene insertion into the plastid genome alters expression of adjacent native chloroplast genes at the transcriptional and translational levels.

In plant biotechnology and basic research, chloroplasts have been used as chassis for the expression of various transgenes. However, potential unintended side effects of transgene insertion and high-level transgene expression on the expression of native chloroplast genes are often ignored and have not been studied comprehensively. Here, we examined expression of the chloroplast genome at both the transcriptional and translational levels in five transplastomic tobacco (Nicotiana tabacum) lines carrying the identical aadA resistance marker cassette in diverse genomic positions. Although none of the lines exhibits a pronounced visible phenotype, the analysis of three lines that contain the aadA insertion in different locations within the petL-petG-psaJ-rpl33-rps18 transcription unit demonstrates that transcriptional read-through from the aadA resistance marker is unavoidable, and regularly causes overexpression of downstream sense-oriented chloroplast genes at the transcriptional and translational levels. Investigation of additional lines that harbour the aadA intergenically and outside of chloroplast transcription units revealed that expression of the resistance marker can also cause antisense effects by interference with transcription/transcript accumulation and/or translation of downstream antisense-oriented genes. In addition, we provide evidence for a previously suggested role of genomically encoded tRNAs in chloroplast transcription termination and/or transcript processing. Together, our data uncover principles of neighbouring effects of chloroplast transgenes and suggest general strategies for the choice of transgene insertion sites and expression elements to minimize unintended consequences of transgene expression on the transcription and translation of native chloroplast genes.

Chloroplast gene expression level is negatively correlated with evolutionary rates and selective pressure while positively with codon usage bias in Ophioglossum vulgatum L.

BACKGROUND: Characterization of the key factors determining gene expression level has been of significant interest. Previous studies on the relationship among evolutionary rates, codon usage bias, and expression level mostly focused on either nuclear genes or unicellular/multicellular organisms but few in chloroplast (cp) genes. Ophioglossum vulgatum is a unique fern and has important scientific and medicinal values. In this study, we sequenced its cp genome and transcriptome to estimate the evolutionary rates (dN and dS), selective pressure (dN/dS), gene expression level, codon usage bias, and their correlations. RESULTS: The correlation coefficients between dN, dS, and dN/dS, and Transcripts Per Million (TPM) average values were -0.278 (P = 0.027 < 0.05), -0.331 (P = 0.008 < 0.05), and -0.311 (P = 0.013 < 0.05), respectively. The codon adaptation index (CAI) and tRNA adaptation index (tAI) were significantly positively correlated with TPM average values (P < 0.05). CONCLUSIONS: Our results indicated that when the gene expression level was higher, the evolutionary rates and selective pressure were lower, but the codon usage bias was stronger. We provided evidence from cp gene data which supported the E-R (E stands for gene expression level and R stands for evolutionary rate) anti-correlation.

Differences in the Chloroplast Genome and Its Regulatory Network among Cathaya argyrophylla Populations from Different Locations in China.

Cathaya argyrophylla Chun et Kuang is a severely endangered, tertiary relict plant unique to China whose high physiological sensitivity to the environment, including photosensitivity, is likely closely related to its endangered status; however, the exact mechanism responsible has remained unknown due to the rarity of the plant and the difficulties involved in performing physiological studies on the molecular level. In this study, the chloroplast genomes of six C. argyrophylla populations sampled from different locations in China were characterized and compared. In addition, a gene regulatory network of the polymorphic chloroplast genes responsible for regulating genes found elsewhere in the plant genome was constructed. The result of the genome characterization and comparison showed that the genome characteristics, the gene composition, and the gene sequence of the chloroplast genes varied by location, and the gene regulatory network showed that the differences in growth location may have led to variations in the protein-coding chloroplast gene via various metabolic processes. These findings provide new insights into the relationship between chloroplasts and the sensitive metabolism of C. argyrophylla and provide additional reference materials for the conservation of this endangered plant.

The Codon Statistics Database: A Database of Codon Usage Bias.

We present the Codon Statistics Database, an online database that contains codon usage statistics for all the species with reference or representative genomes in RefSeq (over 15,000). The user can search for any species and access two sets of tables. One set lists, for each codon, the frequency, the Relative Synonymous Codon Usage, and whether the codon is preferred. Another set of tables lists, for each gene, its GC content, Effective Number of Codons, Codon Adaptation Index, and frequency of optimal codons. Equivalent tables can be accessed for (1) all nuclear genes, (2) nuclear genes encoding ribosomal proteins, (3) mitochondrial genes, and (4) chloroplast genes (if available in the relevant assembly). The user can also search for any taxonomic group (e.g., "primates") and obtain a table comparing all the species in the group. The database is free to access without registration at http://codonstatsdb.unr.edu.

High-resolution map of plastid-encoded RNA polymerase binding patterns demonstrates a major role of transcription in chloroplast gene expression.

Plastids contain their own genomes, which are transcribed by two types of RNA polymerases. One of those enzymes is a bacterial-type, multi-subunit polymerase encoded by the plastid genome. The plastid-encoded RNA polymerase (PEP) is required for efficient expression of genes encoding proteins involved in photosynthesis. Despite the importance of PEP, its DNA binding locations have not been studied on the genome-wide scale at high resolution. We established a highly specific approach to detect the genome-wide pattern of PEP binding to chloroplast DNA using plastid chromatin immunoprecipitation-sequencing (ptChIP-seq). We found that in mature Arabidopsis thaliana chloroplasts, PEP has a complex DNA binding pattern with preferential association at genes encoding rRNA, tRNA, and a subset of photosynthetic proteins. Sigma factors SIG2 and SIG6 strongly impact PEP binding to a subset of tRNA genes and have more moderate effects on PEP binding throughout the rest of the genome. PEP binding is commonly enriched on gene promoters, around transcription start sites. Finally, the levels of PEP binding to DNA are correlated with levels of RNA accumulation, which demonstrates the impact of PEP on chloroplast gene expression. Presented data are available through a publicly available Plastid Genome Visualization Tool (Plavisto) at https://plavisto.mcdb.lsa.umich.edu/.

The role of PAP4/FSD3 and PAP9/FSD2 in heat stress responses of chloroplast genes.

Chloroplasts' mechanisms of adaptation to elevated temperatures are largely determined by the gene expression of the plastid transcription apparatus. Gene disruption of iron-containing superoxide dismutase PAP4/FSD3 and PAP9/FSD2, which are parts of the DNA-RNA polymerase complex of plastids, contributed to a decrease in resistance to oxidative stress caused by the prolonged action of elevated temperatures (5 days, 30 °C). Under heat stress conditions, pap4/fsd3 and pap9/fsd2 mutants showed a decline in chlorophyll content and photosynthesis level, as measured by photosynthetic parameters, and a different amplitude of HSP gene response to heat stress. The expression of nuclear- and plastid-encoded photosynthesis genes and corresponding proteins was strongly inhibited in the mutants as compared with wild-type plants and was further suppressed or displayed no additional changes at 30 °C. NEP-dependent plastid genes, as well as NEP genes RPOTp and RPOTmp, were also downregulated in the mutants by high temperature or remained stable, unlike in wild-type seedlings where these genes were strongly upregulated. The results obtained correspond to the concept of the complex effect of various forms of reactive oxygen species under all types of stresses, including heat stress, and confirm the hypothesis of a new regulatory function in plastid transcription acquired by enzymatic proteins during evolution.

Guanosine tetraphosphate (ppGpp) accumulation inhibits chloroplast gene expression and promotes super grana formation in the moss Physcomitrium (Physcomitrella) patens.

The nucleotides guanosine tetraphosphate and pentaphosphate (or (p)ppGpp) are implicated in the regulation of chloroplast function in plants. (p)ppGpp signalling is best understood in the model vascular plant Arabidopsis thaliana in which it acts to regulate plastid gene expression to influence photosynthesis, plant development and immunity. However, little information is known about the conservation or diversity of (p)ppGpp signalling in other land plants. We studied the function of ppGpp in the moss Physcomitrium (previously Physcomitrella) patens using an inducible system for triggering ppGpp accumulation. We used this approach to investigate the effects of ppGpp on chloroplast function, photosynthesis and growth. We demonstrate that ppGpp accumulation causes a dramatic drop in photosynthetic capacity by inhibiting chloroplast gene expression. This was accompanied by the unexpected reorganisation of the thylakoid system into super grana. Surprisingly, these changes did not affect gametophore growth, suggesting that bryophytes and vascular plants may have different tolerances to defects in photosynthesis. Our findings point to the existence of both highly conserved and more specific targets of (p)ppGpp signalling in the land plants that may reflect different growth strategies.

Beyond RuBisCO: convergent molecular evolution of multiple chloroplast genes in C4 plants.

BACKGROUND: The recurrent evolution of the C4 photosynthetic pathway in angiosperms represents one of the most extraordinary examples of convergent evolution of a complex trait. Comparative genomic analyses have unveiled some of the molecular changes associated with the C4 pathway. For instance, several key enzymes involved in the transition from C3 to C4 photosynthesis have been found to share convergent amino acid replacements along C4 lineages. However, the extent of convergent replacements potentially associated with the emergence of C4 plants remains to be fully assessed. Here, we conducted an organelle-wide analysis to determine if convergent evolution occurred in multiple chloroplast proteins beside the well-known case of the large RuBisCO subunit encoded by the chloroplast gene rbcL. METHODS: Our study was based on the comparative analysis of 43 C4 and 21 C3 grass species belonging to the PACMAD clade, a focal taxonomic group in many investigations of C4 evolution. We first used protein sequences of 67 orthologous chloroplast genes to build an accurate phylogeny of these species. Then, we inferred amino acid replacements along 13 C4 lineages and 9 C3 lineages using reconstructed protein sequences of their reference branches, corresponding to the branches containing the most recent common ancestors of C4-only clades and C3-only clades. Pairwise comparisons between reference branches allowed us to identify both convergent and non-convergent amino acid replacements between C4:C4, C3:C3 and C3:C4 lineages. RESULTS: The reconstructed phylogenetic tree of 64 PACMAD grasses was characterized by strong supports in all nodes used for analyses of convergence. We identified 217 convergent replacements and 201 non-convergent replacements in 45/67 chloroplast proteins in both C4 and C3 reference branches. C4:C4 branches showed higher levels of convergent replacements than C3:C3 and C3:C4 branches. Furthermore, we found that more proteins shared unique convergent replacements in C4 lineages, with both RbcL and RpoC1 (the RNA polymerase beta' subunit 1) showing a significantly higher convergent/non-convergent replacements ratio in C4 branches. Notably, more C4:C4 reference branches showed higher numbers of convergent vs. non-convergent replacements than C3:C3 and C3:C4 branches. Our results suggest that, in the PACMAD clade, C4 grasses experienced higher levels of molecular convergence than C3 species across multiple chloroplast genes. These findings have important implications for our understanding of the evolution of the C4 photosynthesis pathway.

Chloroplast Genes Are Involved in The Male-Sterility of K-Type CMS in Wheat.

The utilization of crop heterosis can greatly improve crop yield. The sterile line is vital for the heterosis utilization of wheat (Triticum aestivum L.). The chloroplast genomes of two sterile lines and one maintainer were sequenced using second-generation high-throughput technology and assembled. The nonsynonymous mutated genes among the three varieties were identified, the expressed difference was further analyzed by qPCR, and finally, the function of the differentially expressed genes was analyzed by the barley stripe mosaic virus-induced gene silencing (BSMV-VIGS) method. A total of 16 genes containing 31 nonsynonymous mutations between K519A and 519B were identified. There were no base mutations in the protein-encoding genes between K519A and YS3038. The chloroplast genomes of 519B and K519A were closely related to the Triticum genus and Aegilops genus, respectively. The gene expression levels of the six selected genes with nonsynonymous mutation sites for K519A compared to 519B were mostly downregulated at the binucleate and trinucleate stages of pollen development. The seed setting rates of atpB-silenced or ndhH-silenced 519B plants by BSMV-VIGS method were significantly reduced. It can be concluded that atpB and the ndhH are likely to be involved in the reproductive transformation of 519B.

Chlorosis seedling lethality 1 encoding a MAP3K protein is essential for chloroplast development in rice.

BACKGROUND: Mitogen-activated protein kinase (MAPK) cascades are conserved signaling modules in eukaryotic organisms and play essential roles in immunity and stress responses. However, the role of MAPKs in chloroplast development remains to be evidently established. RESULTS: In this study, a rice chlorosis seedling lethality 1 (csl1) mutant with a Zhonghua11 (ZH11, japonica) background was isolated. Seedlings of the mutant were characterized by chlorotic leaves and death after the trefoil stage, and chloroplasts were observed to contain accumulated starch granules. Molecular cloning revealed that OsCSL1 encoded a MAPK kinase kinase22 (MKKK22) targeted to the endoplasmic reticulum (ER), and functional complementation of OsCSL1 was found to restore the normal phenotype in csl1 plants. The CRISPR/Cas9 technology was used for targeted disruption of OsCSL1, and the OsCSL1-Cas9 lines obtained therein exhibited yellow seedlings which phenocopied the csl1 mutant. CSL1/MKKK22 was observed to establish direct interaction with MKK4, and altered expression of MKK1 and MKK4 was detected in the csl1 mutant. Additionally, disruption of OsCSL1 led to reduced expression of chloroplast-associated genes, including chlorophyll biosynthetic genes, plastid-encoded RNA polymerases, nuclear-encoded RNA polymerase, and nuclear-encoded chloroplast genes. CONCLUSIONS: The findings of this study revealed that OsCSL1 played roles in regulating the expression of multiple chloroplast synthesis-related genes, thereby affecting their functions, and leading to wide-ranging defects, including chlorotic seedlings and severely disrupted chloroplasts containing accumulated starch granules.

Context-Dependent Mutation Dynamics, Not Selection, Explains the Codon Usage Bias of Most Angiosperm Chloroplast Genes.

Two competing proposals about the degree to which selection affects codon usage of angiosperm chloroplast genes are examined. The first, based on observations that codon usage does not match expectations under the naïve assumption that base composition will be identical at all neutral sites, is that selection plays a significant role. The second is that codon usage is determined almost solely by mutation bias and drift, with selection influencing only one or two highly expressed genes, in particular psbA. First it is shown that, as a result of an influence of neighboring base composition on mutation dynamics, compositional biases are expected to be widely divergent at different sites in the absence of selection. The observed mutation properties are then used to predict expected neutral codon usage biases and to show that observed deviations from the naïve expectations are in fact expected given the context-dependent mutational dynamics. It is also shown that there is a match between the observed and expected codon usage when context effects are taken into consideration, with psbA being a notable exception. Overall, the data support the model that selection is not a widespread factor affecting the codon usage of angiosperm chloroplast genes and highlight the need to have an accurate model of mutational dynamics.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing.

KEY MESSAGE: Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnPGGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

On the Edge of Dispensability, the Chloroplast ndh Genes.

The polypeptides encoded by the chloroplast ndh genes and some nuclear genes form the thylakoid NADH dehydrogenase (Ndh) complex, homologous to the mitochondrial complex I. Except for Charophyceae (algae related to higher plants) and a few Prasinophyceae, all eukaryotic algae lack ndh genes. Among vascular plants, the ndh genes are absent in epiphytic and in some species scattered among different genera, families, and orders. The recent identification of many plants lacking plastid ndh genes allows comparison on phylogenetic trees and functional investigations of the ndh genes. The ndh genes protect Angiosperms under various terrestrial stresses, maintaining efficient photosynthesis. On the edge of dispensability, ndh genes provide a test for the natural selection of photosynthesis-related genes in evolution. Variable evolutionary environments place Angiosperms without ndh genes at risk of extinction and, probably, most extant ones may have lost ndh genes recently. Therefore, they are evolutionary endpoints in phylogenetic trees. The low number of sequenced plastid DNA and the long lifespan of some Gymnosperms lacking ndh genes challenge models about the role of ndh genes protecting against stress and promoting leaf senescence. Additional DNA sequencing in Gymnosperms and investigations into the molecular mechanisms of their response to stress will provide a unified model of the evolutionary and functional consequences of the lack of ndh genes.

Pseudogenization of the chloroplast threonine (trnT-GGU) gene in the sunflower family (Asteraceae).

The chloroplast genome evolves through the course of evolution. Various types of mutational events are found within the chloroplast genome, including insertions-deletions (InDels), substitutions, inversions, gene rearrangement, and pseudogenization of genes. The pseudogenization of the chloroplast threonine (trnT-GGU) gene was previously reported in Cryptomeria japonica (Cupressaceae), Pelargonium × hortorum (Geraniaceae), and Anaphalis sinica and Leontopodium leiolepis of the tribe Gnaphalieae (Asteroideae, Asteraceae). Here, we performed a broad analysis of the trnT-GGU gene among the species of 13 subfamilies of Asteraceae and found this gene as a pseudogene in core Asteraceae (Gymnarrhenoideae, Cichorioideae, Corymbioideae, and Asteroideae), which was linked to an insertion event within the 5' acceptor stem and is not associated with ecological factors such as habit, habitat, and geographical distribution of the species. The pseudogenization of trnT-GGU was not predicted in codon usage, indicating that the superwobbling phenomenon occurs in core Asteraceae in which a single transfer RNA (trnT-UGU) decodes all four codons of threonine. To the best of our knowledge, this is the first evidence of a complete clade of a plant species using the superwobbling phenomenon for translation.

PPR647 Protein Is Required for Chloroplast RNA Editing, Splicing and Chloroplast Development in Maize.

Chloroplasts play an essential role in plant growth and development. Any factors affecting chloroplast development will lead to abnormal plant growth. Here, we characterized a new maize mutant, albino seedling mutant 81647 (as-81647), which exhibits an entirely albino phenotype in leaves and eventually died before the three-leaf stage. Transmission electron microscopy (TEM) demonstrated that the chloroplast thylakoid membrane was impaired and the granum lamellae significantly decreased in as-81647. Map-based cloning and transgenic analysis confirmed that PPR647 encodes a new chloroplast protein consisting of 11 pentratricopeptide repeat domains. Quantitative real-time PCR (qRT-PCR) assays and transcriptome analysis (RNA-seq) showed that the PPR647 mutation significantly disrupted the expression of PEP-dependent plastid genes. In addition, RNA splicing and RNA editing of multiple chloroplast genes showed severe defects in as-81647. These results indicated that PPR647 is crucial for RNA editing, RNA splicing of chloroplast genes, and plays an essential role in chloroplast development.

A Transcription Factor Regulates Gene Expression in Chloroplasts.

The chloroplast is a semi-autonomous organelle with its own genome. The expression of chloroplast genes depends on both chloroplasts and the nucleus. Although many nucleus-encoded proteins have been shown to localize in chloroplasts and are essential for chloroplast gene expression, it is not clear whether transcription factors can regulate gene expression in chloroplasts. Here we report that the transcription factor NAC102 localizes in both chloroplasts and nucleus in Arabidopsis. Specifically, NAC102 localizes in chloroplast nucleoids. Yeast two-hybrid assay and co-immunoprecipitation assay suggested that NAC102 interacts with chloroplast RNA polymerases. Furthermore, overexpression of NAC102 in chloroplasts leads to reduced chloroplast gene expression and chlorophyll content, indicating that NAC102 functions as a repressor in chloroplasts. Our study not only revealed that transcription factors are new regulators of chloroplast gene expression, but also discovered that transcription factors can function in chloroplasts in addition to the canonical organelle nucleus.

Novel Structural Variation and Evolutionary Characteristics of Chloroplast tRNA in Gossypium Plants.

Cotton is one of the most important fiber and oil crops in the world. Chloroplast genomes harbor their own genetic materials and are considered to be highly conserved. Transfer RNAs (tRNAs) act as "bridges" in protein synthesis by carrying amino acids. Currently, the variation and evolutionary characteristics of tRNAs in the cotton chloroplast genome are poorly understood. Here, we analyzed the structural variation and evolution of chloroplast tRNA (cp tRNA) based on eight diploid and two allotetraploid cotton species. We also investigated the nucleotide evolution of chloroplast genomes in cotton species. We found that cp tRNAs in cotton encoded 36 or 37 tRNAs, and 28 or 29 anti-codon types with lengths ranging from 60 to 93 nucleotides. Cotton chloroplast tRNA sequences possessed specific conservation and, in particular, the Ψ-loop contained the conserved U-U-C-X3-U. The cp tRNAs of Gossypium L. contained introns, and cp tRNAIle contained the anti-codon (C-A-U), which was generally the anti-codon of tRNAMet. The transition and transversion analyses showed that cp tRNAs in cotton species were iso-acceptor specific and had undergone unequal rates of evolution. The intergenic region was more variable than coding regions, and non-synonymous mutations have been fixed in cotton cp genomes. On the other hand, phylogeny analyses indicated that cp tRNAs of cotton were derived from several inferred ancestors with greater gene duplications. This study provides new insights into the structural variation and evolution of chloroplast tRNAs in cotton plants. Our findings could contribute to understanding the detailed characteristics and evolutionary variation of the tRNA family.